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We calculate the present expansion of our Universe endowed with relict colored objects - quarks 
and gluons - that survived hadronization either as isolated islands of quark-giuon "nuggets" or 
spread uniformly in the Universe. In the first scenario, the quark nuggets can play the role of dark 
matter. In the second scenario, we demonstrate that uniform colored objects can play the role of 
dark energy providing the late-time accelerating expansion of the Universe. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.+X, 12.39.Ba 



I. INTRODUCTION 

Over two decades ago ^^ (see also Q) the accelerated 
expansion of the early Universe was derived from a quark 
bag model with theproper equations of state (EoS). It 
was called tepid [l|, or little [3| inflation, in view of its 
moderate scales, compared to the better known earlier 
inflation. However, occurring at a later time (when the 
temperature ~ 200 MeV) and smearing a lot of the ear- 
lier effects, it may have important consequence for the 
observable Universe. 

The derivation was based on a quark-gluon bag EoS 
completing the Friedmann equations. Our Universe was 
cooling down along the "hot" (i.e. quark-gluon) branch 
of the EoS until it reached the point/area of the transition 
to the confined hadron phase. In most of the papers along 
these lines (see [1] and references therein) the process of 
inflation terminates by a phase transition (hadronization) 
to the state of colorless objects. The details of the phase 
transition, depending on unknown confining forces, are 
poorly known and leave much room for speculations. 

In the present paper, we consider the possibility that 
a small fraction of colored objects - quarks and gluons 
- escaped hadronization. They may survive as islands of 
colored particles, called quark-gluon nuggets (for brevity 
sometimes also quark nuggets). This possibility was first 
considered by E. Witten [^] and scrutinized further in 
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[^-H]- If SO J the "hot" quark-gluon phase in the form of 
quark nuggets may affect the present expansion of the 
Universe. Indeed, our investigation shows that nuggets 
can contribute to dark matter provided that their inter- 
action with ordinary matter is weak. 

Another possibility is that a very small (to be spec- 
ified!) fraction of colored objects - quarks and glu- 
ons - survived after the phase transition in the form 
of a perfect fluid uniformly spread within the colorless 
hadronic medium. This picture is physically less moti- 
vated than the nugget model. We suppose that this fluid 
has the same thermodynamical properties as the quark- 
gluon plasma (QGP). Therefore, we call it as a QGP-like 
perfect fluid. Nevertheless, it is of interest to investigate 
cosmological consequences of such assumption. In our 
paper, we demonstrate that such fluid can provide an al- 
ternative (with respect to the cosmological constant) ex- 
planation to the late-time accelerating expansion of the 
Universe. 

The paper is structured as follows. In Sec. II, we 
briefly remind the quark bag equations of state. In Sec. 
Ill and IV, we consider the influence of nuggets and QGP- 
like perfect fluid on the late-time expansion of the Uni- 
verse. The main results are briefly summarized in con- 
cluding Sec. V. 



II. EQUATIONS OF STATE 
IN THE QUARK-GLUON BAG MODEL 

We flrst briefly remind the quark bag equation of state, 
a simple model of quark confinement. For vanishing 
chemical potential, // = 0, it is a system of two equa- 
tions 



III. QUARK NUGGETS 



Pq{T)^A,T^~B 



Ph{T) = AuT^ . 



(1) 



(2) 



The first line corresponds to the "hot" phase of decon- 
fined quarks and gluons, and the second one relates to 
confined particles, i.e. hadrons. A system of strongly 
interacting particles, made of free quarks and gluons, is 
cooling down and meets the " cold" phase transforming in 
colorless hadrons. The coefficients are defined by the de- 
grees of freedom and are equal to: A^ w 1.75, A/j w 0.33, 
B^{Aq- Ah)T^ and T^ w 200 MeV. 

Knowing the pressure, p(T), for /i = 0, one can eas- 
ily calculate the remaining thermodynamical quantities, 
e.g., for the energy density we have 



^in-T%-.. 



(3) 



The above EoS is not unique. There is a number of in- 
teresting modifications [i|, [2, Iml2l ■ First such modifica- 
tion was considered by C. Kallmann Q, who introduced 
a temperature-dependent bag " constant" , namely, by re- 
placing in the first fine of the EoS, Eq. ^, B ^ B{T) ^ 
BT, where B = {Aq — Afi)TJ?. This modification has im- 
mediate consequences, namely, by producing a minimum 
in the "hot" line of the EoS, corresponding to metastable 
deeply supercooled states of the deconfined strongly in- 
teracting matter. Also, it drives inflation of the Universe, 
as shown in [l|, 01 . A detailed discussion of the above EoS 
and their consequences, both for the heavy ion collisions 
and the early Universe, can be found in the review paper 

0-. 

Since the idea of the present paper is that a small frac- 
tion of deconfined quarks and gluons survives to present 
days, we shall be interested in the "hot" branch of the 
bag EoS. As we mentioned above, there is a number of 
different modifications of Eq. ([ij. For our present pur- 
poses, however, two simple representatives will be suffi- 
cient. They are the Kallmann modified model (which we 
call Model I): 



PqiT) ^ AqT" - BT = AiT + AiT\ (4) 

and the original model (Model II) described by Eq. ([1]): 



p, (T) ^ A^T" - B = Ao + AiT\ (5) 



It is worth noting that in these equations, we mea- 
sure temperature in energetic units, i.e. in erg or MeV 
(IMeV w 0.1602 x 10~^ erg). Then, pressure is measured 
in erg^ or MeV ^. 



As we wrote above, there is a possibility that after 
a phase transition from quark gluon plasma (QGP) to 
hadronic matter, a part of QGP was preserved in the form 
of quark gluon nuggets |3-0| . They are isolated " islands" 
of QGP in a sea of a new hadronic phase. Now, we want 
to investigate cosmological consequences of this assump- 
tion. Obviously, for considered models, a cosmological 
scenario strongly depends on thermodynamical proper- 
ties of QGP. We focus on two possible Eqs. dU and ([5]). 
With the help of standard thermodynamical Eq. ([3]) we 
get the expressions for the energy density: 



and 



£ = 3A4T* 



£ = -Ao + SA^r 



(6) 



(7) 



for Model I and Model II, respectively. Eqs. (g]), (P, dS]) 
and ([7]) describe the pressure and energy density inside 
of the nuggets. The total pressure and energy density of 
all nuggets in the Universe can be calculated as follows. 
Let us take, e.g.. Model I with Eq. dH). Then, for total 
pressure of nuggets we get 



P = 



V 



AiT + A^T'^ 



(8) 



where pqi is the pressure of the i-th nugget with the vol- 
ume Vi and V (X a^ is the total volume of the Universe 
(a is the scale factor of the Friedmann- Robertson- Walker 
metric) . We consider the case where all nuggets have the 
same pressure ^ and their volumes are either constant 
or only slightly varying with time. The total volume of 
nuggets J2i "^i is included in the coefficients Ai and ^4 
(i.e. ^1,4 have dimension ^1,4 x cm"^, so, taking into ac- 
count the footnote [Tl Ai is dimensionless and A4 has the 
dimension crg""^). Therefore, 



4i = 4^ = -0.8114 r3. 
Ai A4 



(9) 



Similarly, from Eq. ([5]), for the energy density of all 
nuggets we get: 



S = 



3AaT'^ 



(10) 



The same procedure holds for the Model II. Let us con- 
sider two models separately. 



^ Usually, the dimension of pressure is erg/cm'^ . It is not difficult 



to get the relation IMeV* ~ 2.09 X lO^^erg/cm^. However, to 
transform to the usual units, it is more convenient to redefine the 
coefficients as follows: Ai ^ Ai = Ai/[{Mpic^)^L^i] , i = 0, 1, 4, 
where Mpi a 2.177 X 10~^g is the Planck mass and Lpi « 
1.616 X 10~'^'^cm is the Planck length. 



Model I. 



where the cosmological parameters are 



Here, the pressure and energy density of all nuggets are 
given by the above formulae ^ and (ITU)) , respectively. 
In these formulae, temperature is a function of the scale 
factor a: T = T{a). Let us specify this dependence. 
From the energy conservation equation 



d{£a^) + Pd{a^)=Q 



we can easily get 



T = 



Aa 



(11) 



(12) 



As we mentioned above, we consider the model where 
the coefficients Ai < and A^ > 0. In Eq. ([12]), C > 
is the constant of integration which is defined by the 
temperature Tq and scale factor oq at the present time: 

C = (Ai + ^4^0*)'^' «o - At/' (-O.SlUTi* + T^f ao . 

(13) 
The temperature T tends to the constant value when the 
scale factor approaches infinity: 



T 



Tr. 



Aa 



1/3 



0.9327 Tc for a -^ oo , 



(14) 
and the pressure goes asymptotically to zero: P ^- 0. On 
the other hand, for C = 0, the temperature is constant 
all the time T = Too, and nuggets behave as a matter 
with zero pressure P = 0. 

We consider our Universe starting from the moment 
when we can drop the radiation. It is well known that 
the radiation dominated (RD) stage is much shorter than 
the matter dominated (MD) stage. Hence, the neglect of 
the RD stage does not affect much the estimate of the 
lifetime of the Universe. Starting from the MD stage, the 
first Friedmann equation for our model reads 



M' 



K 



«£+«£r*(^)'+A, 



(15) 



where we take into account the cosmological constant A 
and the (usual + dark) matter with the present value of 
the energy density ejf *. In (fTS)) . 'H — a' /a — {da/dr])/a, 
K — SttGn/c^, Gn is the gravitational constant and 
K = ±1,0 is the spatial curvature. The conformal time 
rj is connected with the synchronous time t: adrj = cdt. 
Taking into account Eqs. (|T0| and (rT2|) . we get for the 
Hubble parameter H = {l/a)da/dt = {c/a?')da/dr] the 
following expression: 



H^ = m 



K?)'+H?)' 



n 



M 



(t)^ 



^A + ^K — 



ooV 



a) 



(16) 



VLai — 



3Hf'° 



Qa^ 



ml 



A, 



n 



K 



-K 



( '^ 



and we introduce the dimensionless parameters 

.„2 \3/4 



(17) 



/3 = - 



cV/' 1 



ao/ A', 



A, 



1/4 V am', 



o-'^o 



7 = 



,2 \ 3/4 



A'J' \alHi 



From the second Friedmann equation 

2n' + n^ + K 



kP + A 



(18) 



(19) 



after some obvious algebra we obtain the deceleration 
parameter 



1 d^a 



am dt^ 






\Hj \2 


K?)"*-(?)1 


[H'ir-rCi)""] 


4/3 


^cir-"A- 





1/3 



At the present time to, Eqs. (ITBl) and (1^ read 

1 = {13 + jf + QM + nA + nK, 



(20) 



(21) 



7 



n, 



-90 - -(/3 + 7)'/'-(/3 + 7)'/'-^ + f^A.(22) 



Additionally, we obtain from (1161) the differential equa- 
tion 



dt 



add 



^{P + 7a3/4)4/3 + ^^^a + ^^^4 + ^Ka^ 



where we introduce the dimensionless quantities 

^ a ~ 
a — — , t = Hat . 
ao 



(23) 



(24) 



Therefore, the age of the Universe to is defined by the 
equality 



tn = 



add 



i ^(/3 + 7a3/4)^/^ ^ ^^j^ j^ ^^^4 ^ ^j^i 



(25) 



Now, we consider the case of the flat space K 
VLk^Q- Then, Eq. ^ reads 

! = (/? + 7)'^' + ^M + ^A . 



(26) 



Eq. p2|) demonstrates that accelerated expansion of the 
Universe at the present time (i.e. —go > 0) can be en- 
sured by the first and the last terms on the right side of 
this equation. It is tempting to explain the acceleration 
only at the expense of the first term, i.e. due to the pres- 
ence of quark nuggets when the cosmological constant is 
absent. However, simple analysis of Eqs. (1^^ and ([25)1 
in the case Vl\ = shows that the acceleration — go > 
is achieved only for /3 < that contradicts our model. 
The inclusion of the negative curvature ^k > does not 
affect this conclusion due to the smallness of ^k- 

Nevertheless, quark nuggets can contribute to the dark 
matter if they weekly interact with usual baryon matter 
and light, or they can explain the problem of missing 
baryons [Ij, [ij] if their interaction with usual matter 
is not negligible. As we have mentioned above, nuggets 
behave as matter either asymptotically when a — ;> 00 or 
for all time in the case C = ^> /3 = 0. In the latter case 
we can exactly restore the ACDM model so long as Eqs. 
(l22t and (l26l) take the usual form for this model: 



1 = n 



ALtotal 



^A 



(27) 



and 



v4/3 , 



^A = l-^qo, (28) 



where VtM,totai = 1 ' + ^m- Let VIm correspond 
to just the visible matter. According to observations, 
Q.^[ « 0.04. Then, we can easily restore the parameters 
of the ACDM model. For example, taking the deceler- 
ation parameter go ~ —0.595, as in the ACDM model 



[151, 116|, we get ^a w 0.73 and 7 « 0.33 -^ 7^/2 w 0.23. 
Therefore, ^Mdotai ~ 0.27. For the age of the Universe, 
we get from (^5)1 (where we should put /3 = 0, ^k = 0) 



in 



1 



in 



i/n 



13.7 X 10 yr. Hence, weekly inter- 
acting quark nuggets may be candidates for dark matter. 



B. Model II 

Quark nuggets for the Model II are defined by the ther- 
modynamical functions ^ and ([7]). Similar to Eqs. ^ 
and ([TU]), the total pressure and energy density of all 
nuggets in the Universe are 



where 



p 


Ao + A^^T^ 


(29) 




a3 


E 


-Ao + iA^T^ 
a3 


(30) 


^n 


- '^" - 0.8114 r,^ 
Ai 


(31) 



In this model, Aq < and ^44 > 0. Taking into account 
the footnote [1] we may conclude that the coefficients 
Aq and Ai have dimensions erg and erg~'^, respectively. 
For the thermodynamical functions (1^^ and ([50)1 . energy 
conservation Eq. (fTTj) gives the following dependence of 
the temperature on the scale factor: 



T 



{C/a) - An 
A. 



1/4 



(32) 



where C > is the constant of integration which is de- 
fined by the temperature Tq and the scale factor oq at 
the present time: 

C^{Aq + AiT^) ao = Ai (-0.8114 T,^ + T^) aq • (33) 

Similar to the previous case, the temperature T tends 
to the constant value when the scale factor approaches 
infinity: 



T 



T^ 



-4c 

Ai 



1/4 



0.9491 Tc for a ^ 00 , 



(34) 

and the pressure goes asymptotically to zero: P — )> 0. On 
the other hand, for C = 0, the temperature is constant 
all the time T = Too, and quark nuggets behave as a 
matter with zero pressure P = 0. 

In this model, the pressure and energy density have 
simple dependence on the scale factor: 






£{a) 



3^-44 



(35) 
(36) 



Formally, such perfect fiuid can be considered as a mix- 
ture of radiation and matter. However, for ordinary ra- 
diation T ~ 1/a. 

From the first Friedmann equation ()15|) , we obtain the 
expression for the Hubble parameter: 



H' 



f,h{ 



m 



^M ( 

V a 



^) +7(^)^ 
a / V a / 



Ua 



"^ny 



(37) 



where the cosmological parameters are defined in p7)) 
and the dimensionless parameters /3 and 7 are 






"0^0 



7 



4Ao 
3 






(38) 



Therefore, the age of the Universe to is defined by the 
equality 



-tn = 



ada 



y^/3 -1-70-1- fluja + r^AO^ -I- i^Kl^ 



(39) 



The second Friedmann equation p9)) results in the decel- 
eration parameter 

1 (fa 






f) {-^o'~hm 



l-^i^ 



Qa 



At the present time to, Eqs. ([37| and pO)) read 



1 = /3 + 7 + f^Af + ^A + f^K 

-go = -/3- 



-7 + riA - -f^M. 



(40) 

(41) 
(42) 



It can be easily seen that similar to the previous model 
we also reproduce the standard ACDM model in the case 
of the flat space ^k = and C = -^ /3 = 0. The only 
difference is that in Eqs. ([TT)) and (P5|) i^M, total — 7+^a/- 
For example, if we take f^A/ w 0.04 and go ~ —0.595, 
then we get Ha ~ 0.73 and 7 w 0.23 -^ ^M^totai ~ 0.27 
as in the ACDM model. For these parameters, the age 



of the Universe is to ~ 1 



in 



ifn 



13.7 X lOVr. 



Hence, we again arrive at the conclusion that weekly inter- 
acting quark nuggets may be candidates for dark matter. 



IV. QGP-LIKE PERFECT FLUID 

Above, we considered a scenario where quark nuggets 
form after the phase transition the "isolated islands" in 
the sea of baryon matter. In this section we suppose a 
less physically motivated model where a part of quark 
gluon plasma survived after the phase transition in the 
form of the homogeneously and isotropically distributed 
perfect fluid. This perfect fluid has the thermodynamical 
functions of the form (g]), dH) or (O, ([71). We do not know 
what part of QGP survived (see however some estimate 
at the very end of this section). So, we demand only that 
the ratio between coefficients Ai was preserved. Maybe, 
it is more correct to speak about some unknown perfect 
fluid with the thermodynamical functions motivated by 
the QGP. Therefore, we call this fluid as a QGP-like per- 
fect fluid. We are going to investigate the cosmological 
consequences of such proposal. 

In general, we can consider the pressure of the form 

p{T) ^Ao + MT + MT^ + MT^ + I4T* , (43) 
which, via Eq. ([3]), results in the energy density 

e(r) = -lo + MT"^ + 2A:iT^ + iMT'^ . (44) 

Eq. ([TT|) leads to the differential equation 



^^ ( 2A2 + 6A3T + 12A4T^ 
« 3 f li + 2A2T + 3I3T2 + 4I4T3 



dT . 



(45) 



The solution of this equation enables to determine the de- 
pendence T — T{a). Unfortunately, there is no solution 
of (|45|) in elementary functions. Therefore, we consider 
two particular models by analogy with the previous sec- 
tion. 



A. Model I 

First, we consider the case Aq = 0, A2 = 0, A3 = 0. As 
we mentioned above, the coefficients Ai and A4 satisfy 
the condition similar to ^: 



^^4i = _0.8114r3. 

A4 M 



(46) 



Therefore, we consider the case Ai < and A4 > 0. 
Following the footnote [U the coefficients Ai and A4 have 
dimensions cm~^ and erg~'^cm~'^, respectively. 
Integrating (|45|). we get 



T = 



C^ - Aiq3 
AA^a^ 



1/3 



e(T) = 3A4 



C^ - Aia^ 
4A4a3 



4/3 



(47) 



where C > is the dimensionless constant of integration 
which is defined by the temperature Tq and the scale 
factor flo at the present time: 



d = (^Ai+AA^T^ 



1/3 



flo 



A4^'' (-0.8114 rj' + 4To^ 



.3\l/3 



a-Q . 



(48) 



The difference between the first equation in vf47} and Eq. 
(fT2|) is due to the prefactor 1/a^ in ^ and dTO]). The 
temperature T tends to the constant value when the scale 
factor approaches infinity 



T 



T„o = 



4A4 



1/3 



0.5876 T^ for a ^ 00 . 



(49) 
It can be easily verified that, in the limit a -^ 00, the 
energy density e -^ (3/4)[Af/(4A4)]i/3 and the pres- 
sure p -j> -(3/4)[Af/(4A4)]i/3^ i g_ i-j^g perfect fluid has 
asymptotically the vacuum-like equation of state p — —e. 
If C = 0, then the perfect fluid has this equation of state 
for all time. 

From the first Friedmann equation ([15]), we get the 
Hubble parameter 



H' 



Hi 



K?)^-^' 



4/3 



nMi"-^ 



sia + h 



■< (?)1 ■ 



(50) 



where 



C^ 



4(a3A4)i/4 \at,Hl 



Aial 



3/4 



3/4 



4(a3A;)i/4 \alHl 

We use Eq. ([50)1 to find the age of the Universe: 
_ 



(51) 



I \Jil3 + ja^f/^ + ilMa + fiA54 + 



(52) 



From the second Friedmann equation p^ . the decelera- 
tion parameter is 



H 



27 



"(?)• 



1/3 



V a 



) +7 
a J 



4/3 



- >- (- 

2 V a 



«A \ (i53) 



At the present time ipi Eqs. ([50| and ([5^ read 



(/3 + 7)^/^ + Qm + f^A + f^ 



-90 



K , (54) 

27 (/? + 7)'^' - (/5 + 7)'^' - ^^M + f^A (55) 



The latter equation indicates that the acceleration of 
the Universe (i.e. —go > 0) in this model can originate 
from the first and the last terms on the right-hand side 
of this equation. Up to now, the nature of the cosmolog- 
ical constant is still unclear. So, we try to explain the 
acceleration without it, i.e. we suppose that JIa = 0. It 
is also well known that our Universe is very flat [1^, lla| ■ 
Hence, we put ^k = 0. 

Let us consider two particular cases. The first one 
corresponds to the choice /3 = 0. In this case, Eqs. ([M)) 
and (1551) take the form 



'90 — — :r^M + ^A.gsp , 



(56) 
(57) 



where 51a, ggp = 7'*''^. These equations reproduce exactly 
the standard ACDM model. Therefore, if we take —go ~ 
0.595 and Vtj^i k, 0.27, then we get ^K,qgp ~ 0.73 and for 
the age of the Universe io ~ 1 =i^ ^o ~ -^o^^ ~ 13.7 x 
lOVr- 

The second case with /3 ^ is a bit more complicated, 
but also more interesting. Because the nature of dark 
matter is unclear and, according to the observations [ISl . 
Il6| . the visible matter has J7m ~ 0.04, we shall take 
this value as a total contribution of matter. The other 
experimental restriction follows from the age of globular 
clusters which is 11-^16 Gyr [17|. Therefore, the Universe 
cannot be younger. Now, we shall demonstrate that in 



our model we can satisfy this limitation if we suppose 
only Hm ~ 0.04 and -go ~ 0.595 (as in the ACDM 
model). 

From Eqs. ^^ and ([55]) (where JIa = ^k = 0), we 
express the parameters /? and 7 via f2jv/ and gg: 



2-3r2 



M 



2go 



2-17 



4(l-r!M)i/4 



7 



M 



2go 



4(l-nM)i/4 



(58) 



Then, for Vtu ~ 0.04 and -go ~ 0.595, we get /3 « 0.174 
and 7 w 0.796. For these values of /3 and 7, the age of 
the Universe is to « 0.892 ^ t w 12.2 Gyr, in agreement 
with the experimental data. Roughly speaking, the pa- 
rameter 7 is responsible for the accelerated expansion of 
the Universe, and the parameter l3 plays the role of dark 
substance. 



B. Model II 

Let us consider now the case Ai = 0, A2 = 0, A3 = 0. 
The coefficients Ao and A4 satisfy the condition similar 
to dSB): 



Ao 



A^ 

^4 



-0.8114 Ti'. 



(59) 



where Ao < and A^ > Q^ and they have dimensions 
erg X cm~'^ and erg~^cm~^, respectively. 



Integrating (j45|) . we get 



T 



C 
a 



(60) 



where C > is the constant of integration which is de- 
fined by the temperature Tq and the scale factor ao at 
the present time: C = QoTq. Therefore, the pressure and 
energy density depend on the scale factor as follows: 



e = -Ao + 3^4 — 



^.4 



p = Ao + Ai[ — 



^.4 



(61) 



Hence, such perfect fluid can be formally considered as 
a mixture of vacuum and radiation. This conclusion is 
also conflrmed by the form of the Friedmann equations 
for this model: 



(62) 



H' = Hi 


w-< 




+ n,,{^)\n, + n.{^y' 


(HoV 


hK^)' 


-i.Me)v.A 




where 




V' 


7 


3 \alHl 



(63) 



(64) 



For the age of the Universe we have 



io 



ada 



y^l3 + 7a4 + flMa + fl^a^ + VIko? 



(65) 



Obviously, the paranieter 7 plays the role of the cosmo- 
logical constant. So, we may omit I^a in above equa- 
tions. According to the observations, we may also put 
^K = because of its smallness. We restore exactly 
the ACDM model, e.g., with the choice /3 = (then, 
7 « 0.73). Therefore, in this model the cosmological 
constant arises due to QGP. Let us estimate the fraction 
of QGP which should remain after the phase transition 
to get the observable acceleration. It is clear that param- 
eters ^0 and ^0 play the role of the vacuum energy den- 
sity before and after the phase transition, respectively: 
£v,in = —At), evjin — —Aq. The initial vacuum energy 
density ev,^n = -Aq « 1.42 T* - 2 x lO^MeV^ {%]. For 
7 « 0.73 from (|64|) we get evjin ~ -Aq ^ 6x IG^^erg x 



cm-3 s^ 3 X lG-35MeV*. Thus, syjin/evA 
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V. CONCLUSIONS 

Our paper was devoted to two great challenges of mod- 
ern cosmology and high energy physics dubbed dark mat- 
ter and dark energy. Up to now, there is no satisfactory 
explanation for both of them. In our paper, we proposed 
a possible solution to these problems. For this purpose, 
we considered the expansion of the present Universe, us- 
ing the "hot" , i.e. the quark-gluon branch of the bag 
EoS. Although we made reference to the role of this type 
of the EoS during the early universe, namely its inflation 
phase, here we postponed possible speculations about the 
continuous evolution of the universe, within the present 
formalism, from it early, quark-gluon stage to the present 



days, admitting only the possible continuity in the exis- 
tence in the present Universe of a small fraction of colored 
objects - quarks and gluons - which escaped hadroniza- 
tion. We considered two different scenarios. 

In the first scenario, we supposed that the colored 
objects survived in the form of isolated islands, called 
quark-gluon nuggets, in a sea of a hadronic phase. In the 
second scenario, we assumed that a very small fraction of 
colored objects - quarks and gluons - survived after the 
phase transition in the form of a perfect fluid, called the 
QGP-likc perfect fluid, uniformly spread within the color- 
less hadronic medium. Obviously, these cosmological sce- 
narios are defined by EoS of quark gluon plasma (QGP). 
We focus on two possible Eqs. (|4]) and ([5]) dubbed Model 
I and Model II, respectively. We have shown that within 
considered scenarios, there are no fundamental differ- 
ences in the obtained conclusions for these models. For 
the nugget-scenario, we have shown that weekly inter- 
acting (with visible matter) quark nuggets can play the 
role of dark matter for both of the models. In the case 
of QGP-like fluid, we have demonstrated that this fluid 
can play the role of dark energy providing the late-time 
accelerating expansion of the Universe for both of the 
models. Moreover, we defined that, to be in agreement 
with observations, only 10"^''' part of the colored objects 
should survive after the phase transition. Therefore, the 
considered scenarios provide new possible ways of solving 
the problems of dark matter and dark energy. 
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